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Abstract. We show that the rate of convergence of solutions of finite- 
difference approximations for uniformly elliptic Bellman's equations is 
of order at least /i^", where h is the mesh size. The equations are 
considered in smooth bounded domains. 



Ph ■ The convergence of and error estimates for monotone and consistent ap- 

■^ I proximations to fully nonlinear, first-order PDEs were established a while 

(~| ' ago by Crandall and Lions J5] and Souganidis [24]. 

The convergence of monotone and consistent approximations for fully 
nonlinear, possibly degenerate second-order PDEs was first proved in Barles 
and Souganidis [4j. In a series of papers Kuo and Trudinger [201 \2T\ [22] 
also looked in great detail at the issues of regularity and existence of such 
1^ I approximations for uniformly elliptic equations. 

ir) • There is also a probability part of the story, which started long before see 

^^ I Kushner [18], Kushner and Dupuis J19j . also see Pragarauskas [23]. 

y^ ' However, in the above cited articles apart from [5l[23], related to the first- 

• . order equations, no rate of convergence was established. One can read more 

^— N I about the past development of the subject in Barles and Jakobsen [3] and 

^ ■ the joint article of Hongjie Dong and the author [7]. We are going to discuss 

only some results concerningsecond- order Bellman's equations, which arise 
in many areas of mathematics such as control theory, differential geometry, 
and mathematical finance (see Fleming and Soner [8] , Krylov [9] ) and which 
are most relevant to the results of the present article. 
^_' The first estimates of the rate of convergence for second-order degener- 

ate Bellman's equations appeared in 1997 (see [H]). For equations with 
constant "coefficients" and arbitrary monotone finite-difference approxima- 
tions it was proved in |11] that the rate of convergence is h^'^ if the error 
in approximating the true operators with finite-difference ones is of order 
h on three times continuously differentiable functions. The order becomes 
/i^" if the error in approximating the true operators with finite-difference 
ones is of order /i^ on four times continuously differentiable functions (see 
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Remark 1.4 in [11], which however contains an arithmetical error albeit eas- 
ily correctable. Also see Theorem 5.1 in [Hj). The main idea of [H] that 
the equation and its finite-difference approximation should play symmetric 
roles is also used in the present article. The proofs in [11] are purely ana- 
lytical (in contrast with what one can read in some papers mentioning |llj ) 
even though sometimes probabilistic interpretation of some statements are 
also given. The next step was done in [12] where the so-called method of 
"shaking the coefficients" was introduced to deal with the case of degenerate 
parabolic Bellman's equations with variable coefficients. The two sided error 
estimates were given: from the one side of order h^''^^ and from the other 
h^'^. Here h (unlike in [11]) was naturally interpreted as the mesh size and 
the approximating operators were assumed to approximate the true opera- 
tor with error of order h on three times continuously differentiable functions. 
Until now it is not known whether or not it is possible to improve 1/21 in 
the general setting of [12] . 

However, what is possible is that one can get better estimates if one 
uses some special approximations, say providing the error of order h^ of 
approximating the main part of the true operator on four times continuously 
differentiable functions. This was already mentioned in Remark 1.4 of [11] 
and used by Barles and Jakobsen in [1] to extend the results in [11] to 
equations with variable lower-order "coefficients" . 

One can also consider special finite-difference approximations, for in- 
stance, only containing pure second-order differences in place of second-order 
derivatives, when this /i^ approximating error is automatic. In such cases 
the optimal rate h^'"^ was obtained in the joint work of Hongjie Dong and 
the author [7] for parabolic Bellman's equations with Lipschitz coefficients 
in domains. Both ideas of symmetry and "shaking the coefficients" is used 
in [7] as well as in [6]. In the paper by Hongjie Dong and the author [6] 
we consider among other things weakly nondegenerate Bellman's equations 
with constant "coefficients" in the whole space and obtain the rate of con- 
vergence h, where h is the mesh size. It may be tempting to say that this 
result is an improvement of earlier results, however it is just a better rate 
under different conditions. 

It is worth noting that the set of equations satisfying the conditions in 
[7] is smaller than the one in the papers by Barles and Jakobsen [21 [3], the 
results of which obtained by using the theory of viscosity solutions guarantee 
the rate h^'^. However, in the examples given in [21 [3] of applications of 
the general scheme, for us to get the rate h^'"^, we (only) need to add the 
requirement that the coefficients be twice differentiable (see [H]) and in 
[21 [3] they are only assumed to be once differentiable and still the rate h^'^ 
is guaranteed. One more point to be noted is that in [3] parabolic equations 
are considered with various types of approximation such as Crank-Nicholson 
and splitting-up schemes related to the time derivative. 
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In the present paper we add two more restrictions on the equations from 
[7j: a) we require the coefficients to be in C^'^(M ), b) we require the equa- 
tion to be uniformly nondegenerate. In this case we obtain the rate of 
convergence h?'^, which was announced previously in [13] for equations in 
the whole space. This time "shaking" is not needed as we explain in Remark 
15.11 It would be very interesting to find a way to derive the results of the 
present paper by using methods from [H |2l [3] or other methods based on 
the theory of viscosity solutions. 

1. Main results 

Let A be a set and let 

at = al{x), bt = btix), c" = c"(x), /"(x) 

be real-valued functions of (a, x) defined on AxM.'^ for k = ±1, ..., itdi. We 
assume that, for each x = {xi,...,Xd) G M'', these functions are bounded 
with respect to a £ A. Also let some vectors 

k = ihk, ••., Idk) eB:={x: \x\ < 1} C M'^ 

be defined for k = ±1, ..., itdi. This somewhat unusual range of k turns 
out to be convenient as we explain in Remark 11.11 Consider the following 
Bellman's equation arising, for instance, in the theory of controlled diffusion 
processes (see, for instance, [8], [9]): 

sup [L''v{x) - c"(x)t;(x) + /"(x)] = 0, (1.1) 

aeA 

where 

d d 

L"v{x) = Y, a%{x)Dijv{x)+Yb%x)D,v{x), (1.2) 

A = ^, Ai = A A, 

di d\ 

a%{x) = J2 ^ti^)hkhk, b%x) = J2 bt{x)kk. (1.3) 

|fc|=i |fc|=i 

As follows from the title we will be dealing with uniformly elliptic operators 
and it is well-known that for any uniformly elliptic operator of type ()1.2p 
there exist constant vectors l^ (independent of a) such that representation 
(jl.3p holds. In addition, the regularity properties of afj are inherited by 
a^ (see, for instance, Theorem 11.11 below). This is also known for many 
degenerate elliptic operators (see, for instance, [H]). 

In what follows we adopt the summation convention over all "reasonable" 
values of repeated indices. Observe that 

a%D,,v = a^Dlv, b'^D.v = 6^A,«, 



4 N.V. KRYLOV 

where D^ v and Di^v are the second and the first derivatives of v in the 
direction of Ik, that is 

Df^v = likljkDijV, Di^v = kkDiV. 

Therefore equation (ll.ip is rewritten as 

sup [at{x)Dlv{x) + ht{x)Di^v{x) - c^{x)v{x) + T (x)] = 0, (1.4) 

where, naturally, the summations with respect to k are performed inside the 
sup sign. We approximate solutions of ()1.4p by solutions of finite-difference 
equations obtained after replacing D^ v and Di^v with second- and first- 
order differences, respectively, taken in the direction of l^- 
For any x, ^ G M , /i > 0, and function (j) on M introduce 

Th,mx) =4>[x + h^), dh,ii = — \ ' ^h,i = — ^ j^ ^• 

When ^ is one of the l^s we use the notation 

h,k = ^h,ik^ ^h,k = ^h,ik^ k = ±l,...,±di, 
in which the finite difference approximation of ()1.4p is the following 

sup[aUx)Ah,kvi^) + bt{x)6h,kvi^) " C^^x^x) + T (x)] = 0. (1.5) 

Assumption 1.1. We are given a function g on M.'^ and two constants 
6,K £ (0, go) such that for all a G A and k on M we have 

and for cp = g,a^,b'^,c'^, /", a & A, and j,k € {±1, ..., itdi} we have that 
0gCi.i(Mrf) and 

(7l,l('^l^d^ < K. 



Assumption 1.2. (i) We have 4 = — Lfc, o^ = o°^, A; = ±1, ..., itdi. 
(ii) There exists an integer 1 < lio ^ "^i such that for 

A := {/^,A; = ±l,...,±di}, C := {l±i, ...,l±do} 

we have G £ and 

/: + £ D A D {r + r' : r, i" g £, /'/ ^"}- 

(iii) The coordinates of Ik are rational numbers and Span A = M.'^. 

To justify Assumptions 1 1 . Il and [L2] we remind the reader Theorem 3.1 of 
[15] in which Z = {0, ±1, ±2, ...} and 5^^, 6i > 0, is the set of symmetric 
d X d-matrices a such that for any ^ G M'^ 

<5iiei'<«,o<^r'iei' 

((•, •) stands for the scalar product). 
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Theorem 1.1. There exists a set {li, ..., Z„} C Z such that for any its exten- 
sion {h, ...,lm}, m > n, there exist real-analytic functions Ai(a), ..., Am(a) 
on Ss^ possessing there the following properties: 



m 

E 

fc=i 



^k{a)lklli ^k{a) > 6, VA;, 



where the constant 6 > 0. 



Remark 1.1. Theorem 11.11 implies that any uniformly nondegenerate equa- 
tion of type (jl.ip can be written as (|1.4p with the coefficients satisfying 
Assumption 11.11 as long as the coefficients in (jl.ip are bounded and twice 
continuously differentiable with C^-norm controlled by a constant indepen- 
dent of a and c is uniformly bounded away from zero. If we take {li, ■■.,ln} 
from Theorem 11.11 and define C = {0, Hi, ..., ±Z„} and 

A = {1' + 1" -.I'J" eC,l' ^l"}, 

then Assumption ll.2( ii) will be obviously satisfied and owing to Theorem 
1.11 Assumption 11.11 will still be preserved. Assumption ll.2( i) is of a tech- 
nical nature and easily satisfied just by redefining a^ if necessary which is 
possible since the above A is symmetric and Df = D^^. One could exclude 
Assumption ll.2( i) on the expense of more complicated formulation of As- 
sumption [L2^ii), which, by the way, is needed in order to apply the results 
of |17j about interior estimates of second-order differences of approximate 
solutions. It is also worth saying that by assumption our l^ have lengths 
< 1. Therefore one should normalize those ^^'s from Theorem 11.11 which 
are different from zero, absorbing the extra factors in Afc(a). 

Also observe that Assumption 11.21 requires that Ik = for some k. In- 
cluding the zero vector in A turns out to be convenient from technical point 
of view. The coefficient a^ corresponding to this vector can be set to equal, 
say 1, because the corresponding finite-difference operator is just zero. 

Finally, talking about our assumptions we point out that Assumptions 
1.2( ii)(iii) are used in Section [3] when we apply some results of [16] and [17] 
to finite difference equations in domains. 

We suppose that we are given a "0 S C^(]R'^) such that 

n = {x: i){x) > 0} 

is a bounded domain and l-DV'l > 1 on dQ. Equation (|1.4p is considered in 
0, with V subject to the boundary condition v = g on dU. 
Introduce 

^h = {x e n : x + hB cn}, dhn = M.'^\nh. 

Observe that dh^ contains points which are very far from 0,. This turns out 
to be convenient in our constructions. 

Here are our main results in which the above assumptions are supposed 
to hold and 

p{x) = d[st{x,n^), where 0" = M'^ \ 0. 
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Naturally, p = on Q,'^. 

Theorem 1.2. There are constants N G [0,oo) and /ig > such that for 
all h G (0, Hq] the following is true: 

(i) Equation (jl.Sp in 0,^ with boundary condition v = g on 8^0. has a 
unique hounded Borel solution Vh- 

(ii) On W^ 

P~^\vh - sl, \5h,iVh\, {p - ^h)\5h,ih,jVh\ < N (1.6) 

for any i,j and for any x, y G M*^ 

\vhix)-Vh{y)\<N{\x-y\ + h). (1.7) 

Theorem 1.3. There exists a unique v G Cj^^(J7) n C^'^{]R. ) satisfying 
equation (|1.4p in Q, (a.e.) and equal g on 0^. Furthermore, p\D'^v\ < N in 
i7 (a.e.), where N is a constant. 

This theorem is proved in exactly the same way in which Theorem 8.7 
of [15] is proved. On this way one uses Theorem 11.21 the fact that the 
derivatives of v are weak limits of finite differences of f /i as /i J, (see 
the proof of Theorem 8.7 of [15]), and the fact that there are sufficiently 
many second order derivatives in directions of /j, Ij to conclude from their 
boundedness that the Hessian of v is bounded. 

Remark 1.2. In the theory of fully nonlinear elliptic equations much more 
general results than Theorem 11.31 under much weaker conditions are known. 
For instance, it is proved in [25] that if /3 G (0, 1) is sufficiently small and 
o",6",c",/" G C"{Q) with C^(0)-norms bounded by a constant indepen- 
dent of a, then v G Cf+^(Q) n C^'^{Q). 

The results in [25] and other classical texts on the theory of fully non- 
linear elliptic equations are obtained on the basis of very deep facts, using 
very sophisticated and beautiful techniques, and require a series of long ar- 
guments in the end of which the reader learns a lot of various facts from 
the theory of PDEs and functional analysis. In contrast, our Theorem 11.31 
is obtained on the sole basis of the discrete maximum principle combined 
with elementary albeit quite long computations (see [16] and [T7]) involving 
discrete versions of Bernstein's method. 

Theorem 1.4. There are constants N G [0,oo),/io > such that for all 
h G (0, ho] we have \vh - v\ < Nh"^/^ on R*^. 

Example 1.1. Consider the following uniformly nondegenerate analog of the 
Monge- Ampere equation 

det {-Di,v - ^^5'^Av) = if+y, (1.8) 

where 7 > is a parameter and / G C^'^(M ). It is well known (see, for 
instance, [10]) that equation (11. 8p supplied with the requirement that the 
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matrix {DijV + 'y'^6^^Av) be negative definite is equivalent to the following 
single Bellman's equation: 

sup [{aij + -f^5'^)DijV + {det^/'^a)fd] = 0. (1.9) 

a=a*>0, 
trace a=i 

The above theory is applicable to (jl.Op and shows that we can approximate 
the solutions of (jl.Sp with a C^'^ boundary condition satisfying {DijV + 
7^(5*-' Av) < with solutions of corresponding finite-difference equations. 
It is to be noted however that the number di related to the number of 
directions needed to write (jl.9p in the form (|1.4p will depend on 7 and will 
go to infinity as 7 — )■ 0. 

Also observe that one can prove that if fi is a strictly convex domain, and 
v"^ are solutions of ()1.9p with zero boundary condition, then |t>'^ — t>^| < A^7, 
where v^ is a (generalized or viscosity) solution of (jl.Sp with 7 = 0. 

The rest of the article is organized as follows. In Section [2] we prove 
(jl.7p assuming (]1.6p for equations more general than Bellman's equations. 
Section [3] contains two results from |16| and [T7] needed to prove ()1.6p in 
Sectional In the final rather long Section [5] we present the proof of Theorem 

We use N to denote various constants which may change from one ap- 
pearance to another. Sometimes we specify what they are depending on or 
independent of. However sometimes we do not do that. In these situations 
it is understood that they are independent of anything which is allowed to 
change like x,h,e.... 

2. On the Lipschitz continuity of Vh 

The setting and notation in this section are different from the ones in 
Section [TJ Let 

A:={ik;k = l,...,di} 

be a symmetric subset of B = {x : \x\ < 1}. 

Let H{p, X, u) be a real-valued function given for 

Fix two constants 6 € {0,1], K £ [0,oo). 

Assumption 2.1. The function H{p,x,u) is locally Lipschitz continuous 
with respect to {p,u). Furthermore, at all points of differentiability of H 
with respect to {p, u) we have 

6<H^'j<K, j = l,...,di, H^'^<-6, 
\H,,,\<K, j = 0,...,di, \Hp{p,x,u)\<K{l + \u\). 
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Introduce 

H{x,u) = H{0,x,u). 

The following lemma is an easy consequence of the mean value theorem 
in the form of Hadamard. 

Lemma 2.1. For any p E R, x S M , and u,v £ M. '^'^^ there exist numbers 
ai,..., ad^, bi,..., bdi, c, and f such that 

H{p, X, u) — H{x, v) = H{p, X, u) — H{p, x,v) + f 
= ai[u'l - v'l) + b,{u'i - v[) + c{u'q - v'q) + /, 
5<ai<K, -K<c<-6, \bi\ < K, \f\ < K\p\{l + \v\) (2.1) 
for all i £ {1, ..., di}. 

For any function v and h > define 

Hh[v]{x) = H{x,v{x),5hv{x),Ahv{x)), 

where 

ShV = {6h,iv, ..., 6h,div), Ahv = {Ah,iv, ..., Ah4-^v), 

Theorem 2.2. Let D be a domain in M and let h G (0,5/ A']. Then there 
exists a constant Nq depending only on 6,K, and di such that for any two 
functions v' and v" given on W^ we have that in D^ 

No{v' - v") + h'{H^,[v'] - H,,[v"]) < {No - Sh') sup(7;' - v")+, (2.2) 

sup(u' - v")+ < 5~^ sup{Hh[v"] - Hh[v']), + sup(i;' - v")+. (2.3) 

Furthermore, if H{x, 0) is bounded and we are given a bounded function 
g on W^, then the equation 

Hh[v] = (2.4) 

in Dfi with boundary condition v = g on d^D (in case d^D ^ f/i) has a 
unique bounded solution. This solution is Borel measurable in x if g and 
H{x,u) are Borel for each u. 

Proof. As is usual for monotone finite-difference equations in the proof of 
solvability of (12. 4p we rely on the method called Jacobi iteration in |21] and 
the Banach fixed point theorem in the space *B of bounded functions on M 
provided with the sup norm. We will follow an argument from [15] where 
the situation is somewhat different. 

First we deal with (12. 2p . By Lemma 12. 1^ for w = v' — v" we have on M 
that 

di d\ 

h\Hh[v']{x) - Hh[v"]{x)) = {ch" -hJ2bj - 2^a,)u;(x) 
+hbjw{x + h£j) + aj[w{x + hij) + v{x — hij)]. 
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Since aj, bj, — c < K, and h £ (0, 6/K], there is a constant Nq depending 
only on 5, K, and di such that 

di di 

No + ch^ - h'^bj - 2Y^ ttj > 0. 

Furthermore, 

/i6ji(j(x+/i£j)+aj[it;(x+Mj)+t/;(x— /i^j)] = {hbj+aj)w{x+h£j)+ajw{x—Mj) 

and ah the coefficients on the right are nonnegative in hght of the fact that 
h £ (0, S/K] and aj > 6. It follows that in Dh the left-hand side of (12. 2p is 
less than 

di di 

(A^o + ch — h y^ bj — 2 y^ ttj) sup i(j+ 

di di 

+ y {hbj + ttj) sup w+ + y aj sup Wj^ 

= (iVo + c/i^) sup(i;' - v")+ < {No - 5h^) sup(i;' - v")+, 

D D 

which proves (j2.2p . 

While proving (j2.3p we may assume that 



In that case 



sup(u' — w")+ > sup(u' — v") 



sup(f' — «")+ = sup(t)' — v")j^ 



and ()2.2p implies 

No{v' - v")+ < {No - 5h^) sup(«' - t>")+ + h'{Hh[v"] - Hh[v'\) 

After that it only remains to take the sups of both sides over D^^ and collect 
like terms. 

Now we prove the second assertion of the theorem. Observe that, as is 
easy to see, equation ()2.4p in Dh with the boundary condition v = g on dhD 
is equivalent to the following single equation: 

v{x) = [N^^h^Hh[v]{x) +v{x)\Id^{x) + g{x)IdM^)- (2-5) 

Then introduce an operator T^ : v ^ T^v by 

Thv{x) := [N^^h^Hh[v]{x) + v{x)\Id^{x) + g{x)lQ^D{x). 

Since H\^] is bounded by assumption, TO is bounded. By (12.2p we have that 
ThV = ThV—ThO+ThO is bounded if u is bounded, so that T^ is an operator in 
*B, and moreover Th is a contraction operator in *B with contraction constant 
less than 1 — Nq 5h?' < 1. By the Banach fixed point theorem, equation 
(j2.5p has a unique bounded solution. This solution can be obtained as the 
limit of (r/i)"0 as n — 7- oo. Furthermore, T^ maps Borel functions into Borel 
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measurable ones if H{x, u) is Borel with respect to x. Therefore, under this 
condition, given that g is Borel, all functions T^O are Borel measurable, so 
that the solution v is also Borel measurable. The theorem is proved. 

By taking first v' = v, v" = and then v' = 0, v" = v we obtain from 
([23]) the following. 



Corollary 2.3. If v is a solution of (12.41) with boundary condition v = g 
on dhD, then 

supf+ < 5^^ sup (i?/i[0])+ + sup5+, supf_ < 5""*^ sup (i?^/i[0])_ + sup 5-. 

Rd Dh dhD Rd Dh dhD 

We improve these estimates in Lemma [2. 5 1 for more restricted range of h. 
We need a version of Lemma 8.5 of [15]. Let il. be the set introduced in 
Section [TJ 

Lemma 2.4. Assume that Span (A) = M . Then there exist Hq € (0, oo) 
and a nonnegative function ^ G C^(Jl) such that '^ / p and p/^ are hounded 
in n, for h € (0, /iq] on 0^ we have 

m 

ajA;,,^-^ + kY,\K3^\ < -1, (2.6) 

i=i 

and such that there exist constants M,fj, > 1 for which the function <I> : = 
M$ln(M/^) satisfies 

m 

ajAh,j^ + kY,\K3^\ < -P~' (2.7) 

i=i 
in ilfj,h, whenever 5 < Oj < K . 

Below by Hq we mean a constant in (0, 5/K] for which the statement of 
Lemma 12.41 is true. 



Lemma 2.5. Assume that Span (A) = M'^. Let D d Q,, h ^ (0, /iq] and 
let V satisfy (|2.4p in Dh and v = g in dhD, where g S C^'^{M.'^). Then 
\v — g\ < Np on M. , where the constant N is independent of h. 

Proof. It suffices to prove that |i; — ^rj < N"^ on W^ if we continue ^ 
outside Q as zero. In order to do this we observe that by Lemma |2. II with 
p = at each point of D^ 

Hh[g + N^] = Hh[g + N^] - Hh[0] + Hh[0] 

= NiojAh^i"^ + bjdhj'^ + c^) + ajAh,ig + bjdhjg + eg + H{-, 0), 

where aj,bj,c are some numbers satisfying (|2.ip . Owing to (12. 6j) the last 
expression is negative if we take N large enough. This implies that v — g < 
N^ by the maximum principle (that is by (12. 3p ) and by the fact that v = g 
and ^ > in dhD. The estimate v — g > —N'^ is proved similarly. The 
lemma is proved. 
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In the rest of the section we assume that 

Span (A) = M'^, 

fix a constant /x > 1, and take a constant v G [0,oo). Define 

K = fiV {2u) 

and let D be an open subset of ri^/i- For each p G M and h > introduce 
an operator H^ acting on functions v given on M"^ by the formula 

HhMix) = H{p,x,v{x),6hv{x),Ahv{x)). 

Theorem 2.6. Assume that for each h G (0, Hq] and each p we are given a 
bounded function v^ on M which satisfies the equation 

HM] = (2.8) 

in D^. Introduce 

Vh = vl 
and assume that there is a constant Nq G [0, oo) such that for h G (0, Hq] 
and X G D^ we have 

\vhix)\, \Sh,iVhix)\, {p{x) - iyh)\Ah,iVh{x)\ < No 
for all i. Then for h G (0, /iq] 

Kix) - Vh{x)\ < N{\p\ + sup \vl - Vh\) (2.9) 

dhD 

for any x G M'^ and p G M, where the constant N is independent of x,p, 
and h. 

Proof. Take a constant A'^i to be specified later and introduce 

Wh = '^h + Nip<^. 

By Lemma l2. II for each x 

HI[wh]{x) - H{vh]{x) 

= NlpajAh,^^x) + Nipbj5h,j<^{x) + Nipc^{x) + /, (2.10) 

where aj,bj,c are some numbers satisfying (|2.ip and / is such that 

I/I < K\p\{l + Y,i\^h,Mx)\ + \Sh,Mx)\) + |^;,(x)|). 
j 
Observe that the second term on the left in ()2.10p vanishes in Df^. Further- 
more, < (p — i//i)~^ < 2p~^ if p > 2i//i, so that, owing to the fact that 
K>2u and Dj^ C rj^h, in D/j we obtain 

I/I < N2\p\p-\ 

where N2 is independent of h and p. Finally, we set A^i = A''2, take into 
account (12. 7p and the fact that D^ C ^^h (since k > /i), and we conclude 
from ([2111 that 

Hl[wh] < (2.11) 
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in Dh if p > 0. Upon comparing this with ()2.8p and using (12 .30 we get that 
for any x 

vf^{x) < Wh{x) + sup(?;jj - Wh)+ 

dhD 

<Vh{x)+Np + snv{vl-Vh)+, (2.12) 

dhD 

where A^ is independent oi x,h and p. 

If p < the inequahty in (I2.1ip is reversed and one gets 

Vh- Np<Wh<vl + supi^Wh - vl)+ <vl + sup{vh - vl)+. 

dhD dhD 

By combining this with ()2.12p we come to ()2.9p and the theorem is proved. 

Corollary 2.7. Assume that H{x,u) is locally Lipschitz continuous with 
respect to (x, u) and at all points of its differentiability with respect to {x, u) 

\H^,{x,u)\ <K{l + \u\) 

for all i. Let H[x,0) be bounded on M . For h G {0,5/K] denote by Vh a 
unique bounded solution of ()2.4p in Qh with boundary condition v = g on 
dh^ (see TheoremlKW, where g G C^'^M''). 

Finally, assume that there is a constant Nq £ [0, oo) such that for h £ 
(0, ho] in O/i we have 

I Sh,iVh {x)\, {p{x) - iyh)\ Ah^iVh (x) I < No 

for all i. 

Then for all h G (0, ho] and x,y £ W^ 

\vh{x)-Vhiy)\<N{\x-y\+h), (2.13) 

where the constant N is independent of x,y, and h. 

Proof. If |x — y| > h one can spht the straight segment between x and 
y into adjacent pieces of length h combined with a remaining one of length 
less than h. This shows that we need only prove (j2.13p for |x — y| < h. 

Then fix a unit vector / G W^ and for p G M redefine H if necessary by 
setting 

H{p,x,u) = H{x + pl,u) if \p\ < h, 

H{p,x,u) = H{x + lsignp,u) if \p\ > h, 
that is 

H{p, X, u) = H{x + l(f){p),u), 
where (j){p) = {—h) V p Ah. Observe that the function 

vl{x) := Vh{x + l(j){p)) 

satisfies (12. Sp (with new H) in il.2h ^ ^(k+i)/h where the inclusion follows 
from the fact that k > p > 1. By Lemma [2. 51 there is a constant A'^i G (0, oo) 
such that \vh — g\ < Nip for h G (0, ho]- 
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Now set D = ^Kh- Then Dh = ^(K+i)h a-nd we infer from Theorem 
that for IpI < h 

\vh{x + lp) -Vh{x)\ < N{\p\+ sup \vhi- + Ip) - Vh\) 

^'(K+l)h 

<N{\p\+ sup \vh-g\+ sup \vh-g\+ sup \g{- + lp) - g\). (2.14) 

QC QC OC 

Since |f/i — 5I < -/Vi/O, bfc — 5I < A^i(k + 2)/i outside r2(re+2)h) and this along 
with (|2.14p and the arbitrariness of/ proves (|2.13p for |a;— y| <h< ho, which 
finishes proving the corollary. 

3. Two RESULTS FROM [16] AND [T7] 

We suppose that the assumptions in Section [1] are satisfied and take "cut- 
off" functions 

r?GC2(R'^), |7y|<l, C = ^'• 
Fix an /i G (0, S/{2K)] and set 

Ah,i = hA, Afe,„+i = Afe,n + hA, n > 1, A^^oo = |J A^^^ , 

n 

Define 

Q° = {x£ A/,,00 : 2; + 3hB C ft} = A^^oo n n^h, 
Q = {x + Ah,2-x£Q°}, 6Q = Q\Q°. 

Observe that Q C il-h and Q is a finite set. The latter is due to Assumption 
ll.2( iii) and follows from the fact that the number of points with integral 
coordinates lying in a bounded domain is finite combined with the fact that 
there is a number M such that the coordinates of all points in MAi^oo are 
integers. 

The following is a specification of Theorem 1.2 of [16] in the present 
setting. 



Theorem 3.1. Let v satisfy (II. 5p in Q. Then there is a constant Ni > I 
depending only on di,K, and 6 such that, for any v satisfying 

V > iVi(sup \D'^r]\ + sup \Dr]\'^ + 1), 

we have 

1 A'"i A''i 
iQa.x(\5h,kv\ < Nl{^/l' H — ) max \v\ H — ^ H \- Ni inax(\5h,kv\- (3.1) 

k,Q ' V Q yv V k,SQ 

Remark 3.1. In J16j a more general statement than Theorem 13.11 is proved 
under the assumption that a",6",c°,/'^ are in C'^'^(M'^) rather than in 
C^'^(R'^). Also 6Q in (j3.ip is replaced with a "thinner" set and Assumption 
ll.2( ii) is not used. 

Now comes a version of Theorem 1.1 of 1171. 
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Theorem 3.2. Let v satisfy (jl.Sp in Q. Then there exists a constants N2 
depending only on K, di , and 5 such that if a constant 1^ > satisfies 

V > N2{sup \D'^ri\ + sup \Dr]\'^ + 1), 

then 

T^axC\6h,i6h,jv\ < may: C\6h,iSh,jv\ 

Q,t,j SQ,t,j 

+N2{u^^^ + sup|i:>r/|)max|(5/,,jw| + N2{1 + max|i;|). 

IRd Q,i ' Q 

4. Proof of Theorem 11.21 
(i) For x G M and u = {u',u"), where 

u' = {u'_a^, ...,u'„^,Uq,u[, ...,u'aj, u" = ('u±i,...,'u±^J, 
introduce 

H{x,u) = supK(x)< + 6^(x)< - c^[x)u'^ + f^{x)]. 

By using the inequality 

supF" - supG" < sup(F° - G") (4.1) 

0.0.0. 

we easily conclude that H{x,u) is locally Lipschitz continuous with respect 
to {x,u). As such it is almost everywhere differentiable. Its derivatives are 
limits of finite differences and by using ()4.ip again we see that at all points 
of differentiability of H we have 

\H^.\,\H^:,\<K, j = ±l,...,±di,i = 0,±l,...,±di, H^,<-5. 

Furthermore, 

supF" - supG" > inf(F" - G"), 

a a '^ 

which implies that H^n > 5. 

We have checked that H (independent of p) satisfies Assumption 12.11 
Now assertion (i) follows from Theorem 12.21 and the fact that H{x,u) is a 
(Lipschitz) continuous in x, the latter being again a consequence of ()4.ip 
and Assumption 11.11 This proves assertion (i) and combined with Lemma 
12.51 shows that there is a constant A^ G (0, 00) such that \vh — g\ < Np for 
h £ (0, /io], which is part assertion (ii). 

(ii) In light of (14. ip at all points of differentiability of H{x, u) with respect 
to {x,u) 

\H:,A^,u)\ <N{K,di){l + \u\) 

for all i. Hence, owing to Corollarv 12.71 to prove assertion (ii), it suffices to 
prove (|1.6p . in which the first estimate is obtained above. 
Note that if x E ^h,oo and x Q°, then p{x) < 3/i, and 

\vh{x + hlk) - g{x + hlk)\,\vh{x) - g{x)\ < Nh 
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for any k. Hence for x £ Ah^^o \ Q° we have 

\Sh,kVh{x)\ < h^'^lvhix + hlk) - Vh{x)\ 

< h^Hlvhix + hlk) - g{x + hk)\ + \vh{x) - g{x)\) + N < N, 
where N is independent of h and x. In short 

\Sh,kVh\ < N (4.2) 

on A/i^oo \ Q°, where A^ is independent of h. 

By Theorem 13.11 with t] = 1 for all sufficiently small h 

max \6h,kVh\ < A^(l + max \vh\ + max \5h,kVh\), 

k,Q Rd k,5Q 

where A^ is independent of h. Since \vh — g\ < Np and ()4.2p holds on Afi^^o \ 
Q°, we conclude that ()4.2p holds on A/^oo) provided that h is sufficiently 
small, where A^ is independent of h. Actually, (14. 2p holds on M just because 
any x S R can be placed into an appropriate shift of A/^^oo with the shift 
not affecting any of the above constants A^. 

Thus, we estimated the second quantity in (II. Gh . 

To estimate the last one we use Theorem 13.21 Once again it suffices to 
concentrate on points in Ah,oo- Introduce 2po = max p. It is a standard fact 
that for any r G (0,po] there exists an r)^'^> € Cq°{Q) such that 

rj^^' = 1 on Q2r, f] = outside O^, 

|r/W| < 1, |Z)r/W| < N/r, jD^^^'^)] < N/r"^, 

where the constant A^ is independent of r. Furthermore, since Q" = A^j o^ ^ 
0,2fi and Q C r^ft it holds that SQ C ^\ O,^^- By taking r > 3h, r < 1, 
and applying Theorem 13.21 with t] = t]^^') (and u = N/r"^, where A^ is an 
appropriate constant independent of r and h), for sufficiently small h > 
and any i,j we obtain on A/^oo that 

{v^'Hx)f\6h,^6H,M^)\<N/r, (4.3) 

where A^ is independent of x, r, and h. 

Now we use the arbitrariness of x and r. Take an x £ A^^oo such that 
p{x) > 6h and take r = p{x)/2 (> 3/i). Then r]^'^'{x) = 1 and (|4.3p yields 
p{x)\5h,i5h,jVh{x)\ < N implying that 

{p{x) - 6h)\6h,iSh,M^)\ < N (4.4) 

whenever x G A/,,^oo and p{x) > 6h. However, ()4.4p is obvious if p{x) < 6h. 
Thus ()4.4p holds on Ah^^o and, as it was mentioned, on M . The theorem is 
proved. 
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5. Proof of Theorem 11.41 

We need two additional auxiliary results. Estimate (j5.3p is quite similar 
to (4.6) of J6] the latter being a particular case of the former which occurs 
when N* = (in (4.6) of J6] it is assumed that u € C*''^). 

Lemma 5.1. Let vector I ^ B be on the first basis axis in M . Let e,h > 0, 

T] € C^(M ) be a spherically symmetric function with support in B^ = {x : 
\x\ < e} and let u be a bounded Borel function on M . Define 

B'^ = {z€ R"^-^ : \z\ < e}. 

Assume that for all x,y & [—2h — e,2h + e] x B'^ we have 

\u{x)-u{y)\<Ni\x-y\+N^h, |A/,,;n(x)| < iVa. (5.1) 

Introduce w = u* r]. Then 

\Dfw{0) - Ah,iwm < N2h^\\Df7]\\L, + h\Nie + 7V*/i)||I)f??||L,. (5.2) 

Also for X = 7/, l^l < h, we have 

\dM^)\ < N2\H\l, + h\Nie + Nlh)\\Dfr^\\L,^ (5.3) 

Finally, 

lA^(O) - 6h,iwm < ^2/i||r/||Li + h\Nie + Nlh)\\Dfr]\\L,. (5.4) 

Proof. First for n = 1,2, ..., 6 and x = 7/, |x| < 2/i, we have 



D'^w{x) = u{x- y)D^ri{y) dy = / [u{x - y) - u{x)]Dfr]{y) dy. 

J Be JBe 

Hence, 

\Dfw{x)\ < {N,e + Nlh)\\DYr^\\L,. (5.5) 

Next observe that for smooth functions f{t) of one real variable t we have 

^„(„) , /W- 2/(0) +/(-/.) _ ^ j\, _ |,|,3;<.,(,^),, 

Hence 

\Dfw{Q) - /^h,iw{Q)\ < h^ sup \Dfw{ltih)\. (5.6) 

|ii|<i 

Therefore, to prove ()5.2p it suffices to show that for |ti| < 1 

\Dfw{ltih)\ < N2\\Dfr^\\L,+h\Nie + Nlh)\\D^r]\\L,. (5.7) 

To this end we use (15. 6p which implies that 



\Dfu * r/(0)| < \Ah,iu * r/(0)| + /i^ sup \Dfu * Dfr]{ltih)\. (5.8) 

|ti|<i 

By applying (|5.8p to Dfrj in place of rj and Itih in place of we find 

\Dfw{ltih)\ = \Dfu*Dfr]{ltih)\ < \Ahju* Dfr]{ltih)\ 

+h^ sup \Dfu*Dl7]{l{ti+t2)h)\, 

\t2\<l 
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Here, owing to (I5.ip and the fact that \ti + t2\h + e < 2h + e, 

\Ah,iu*Dfrj{l{ti+t2)h)\ < N2\\Dfri\\L,. 
Furthermore, by (j5.5|) and the fact that |ti + t2\h < 2h 

\Dfu * Dfr]{l{ti + t2)h)\ = \u * Dfr]{l{ti + t2)h)\ 

<{Nie + Nlh)\\Dfrj\\L,. 

This proves ()5.7I) and ()5.2p . 

Since, as above, \Ah^iw{0)\ < N2\\r]\\Li, we obtain (|5.3p for x = directly 
from ()5.6p and ()5.5p . For other values of a; we obtain ()5.3p upon observing 
that the above argument is valid if we replace and Itih in (j5.6p with x and 
rr + /ti/i, respectively. 

To prove (j5.4p observe that in the one-dimensional case 

/(o) = m_M -hj\i- t)f{th) dt. 

It follows that 

\Diw{d) - 5h,iw{<d)\ <h (1 - t)\Dfw{lth)\ dt < h sup \Dfw{lth)\, (5.9) 

./o te[o,i] 

which along with ()5.3p yields (j5.4p . The lemma is proved. 

Remark 5.1. Take a nonnegative spherically symmetric ^ £ C^(M ) whose 
integral is one and whose support is in B and for e > and locally integrable 
functions u on M use the notation 

(^)(x) = / u{x-ey)C{y)dy. 



Denote by Vh a unique bounded solution of (II. 5p . with zero boundary 
data on dhfl. By Theorem 11.21 the function Vh is well defined at least for 
sufficiently small h > 0. 

We want to explain why we iterated ()5.8p and why only once. We will 
apply Lemma l5. II to u = Vh- For simplicity assume that 6 = and notice 
that on Qh for any a £ A we have 

a'iDlv-c''v + r<0. (5.10) 

Next assume that a and c are independent of x as in [6] . Then we can mollify 
all terms in (jS.lOp and obtain that in flh+e we have 

a^A,,4^)-c-^i^) + [r](^)<0. (5.11) 

It is well known that | [/"]('') - /"| < Ne"^ since /" G C^'^ Therefore, (f5lT]l 
implies that 

Now we replace Ah^^v}^ with Df v^ . Ignoring for a moment that our 
estimates of A/j ^ti/i are local in il., we obtain from Lemma |5. II that 

a'^Dlvl^^ - c"?;["^ + /" < iV[e2 + hh^^ + h^e + h)e^^] =: NM{h, e). 
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It follows from Lemma 12.41 that for, perhaps, another constant A*" 

H[vl^^ + N'^M{h,e)] <0 
in 0,h+£ By the maximum principle 

V <vl^^ +NM{h,e)+ sup(i;-w[^^)+. (5.12) 

It seems that the best possible estimate of the last term is that it is less 
than N{h + e) and we will prove this estimate later. Then 

v<v'^^'^ +N[M{h,e) + h + e], 

v<Vh + N[M{h, e) + h + e]+ sup(t;f ^ - Vh)+. (5.13) 

Let us ignore the contribution of the last term in the right-hand side and try 
to make M{h, e) + h + e as mall as possible on the account of arbitrariness 
in choosing e. Observe that this quantity contains e + h?e~''^ which is bigger 
than '^h?/^ , where the constant 7 > is independent of h. Furthermore, 
e + h?e~'^ = 2/1^/3 when e = /i^/^. With this e 

as well and we obtain that v <Vh + Nh?'^ for sufficiently small h > 0. 

This is roughly the way we are going to use (j5.2p . 

Now imagine that we did not enhance the estimate of Dfw and after (|5.6p 
just used (|5.5p to obtain 

\Dfv^^^-A,,ivl:^\<Nh\e + h)e-\ 

Then we would obtain ()5.13p with 

M{h,e) =e^ + h'^{e + h)e-^ 

and M{h, e) + h + e would contain the term e + h?e~'^ whose minimum with 
respect to e is of order h^'"^ and is comparable with its value at e = h}'"^ . 
Then at best we would have that v <Vh + Nh^'"^ . 

On the other hand, we could iterate (j5.8p one more time and obtain that 

\Dfu * Dfr]iliti + t2)h)\ < \Ah,iu * AS(K*i + *2)/i)| 

+/i2 sup \Dfu*Dfri{l{ti+t2 + t3)h)\, 

|t3|<l 

which in our case leads to 

\Dfvl^^ - Af,^ivi'^\ < Nh^[e~^ + h^ie~^ + h\e + h)e~')] 

= N{h^e"^ + /i^e"^ + h^{e + h)e~^). 
This time again M{h, e) + h + e contains e + /i^e~^, which will not lead to 
a better rate than h?'^ . 

It is seen from the above that the fact that |[/"]^'^'* — /"| < Ne^ plays 
no significant role. The estimate K/"]^*"^ — /°| < Ne would do equally well. 
Also in this framework we will be satisfied with estimating the last term in 
the right-hand side of (J5.13P just by iVe, which is quite easy. 
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It is worth noting that the situation with constant in x coefficients a, 6, 
and c is quite different in the whole space (see [6] ) . There no boundary term 
hke the last term in (j5.12p appears and one ends up with v < v^ +NM{h, e). 
With some additional effort (even in our setting of bounded smooth Q and 

variable coefficients, see Lemma [5. 3p one can prove that \vj^ —Vh\ < N{e'^ + 
h). Then 

v<Vh + N[h + e^ + h^e-^ + /i^(e + h)e-^]. (5.14) 

The minimum of e^ + /i^e~^ with respect to e is proportional to its value at 
e = h^''^ and is of order h. Other error terms on the right in (|5.14p are of 
the same or higher order. In this way it is proved in [6] that v < Vh + Nh. 
It is also shown there that in general this estimate is optimal, so that there 
is no need to even consider additional iterations of (15. Sp . 

Finally, we point out that even for the equations in the whole space with 
variable coefficients the error term of order e still appears in the transition 
from (j5.10p to ()5.1ip . The method of "shaking the coefficients" produces 
an error of the same order, which allows us not to use this method on the 
account that we have a good control of the second-order differences of v^ ■ 

Lemma 5.2. There exists a constant N such that for all sufficiently small 
h > 0, for any e > 0, 

(i) In W^ we have \v^^'' - Vh\ < N{e + h); 

(a) In ri4£_|„i6/i for any a ^ A we have 

a^A>P + htDi^v^^ - c^vi\x) + /" 
<N{h + e + h^e-^)p~^ + Nh^e-^{e + h){e^ + h). (5.15) 

Proof. Assertion (i) follows immediately from the fact that \vh{x) — 
Vh{y)\ < N{\x — y\ + h) (see Theorem II. 2p . 

(ii) Fix an a £ A and observe that for x G ^h+e 'we have 



[atAh^kVh + bl6h,kVh - c'^Vh + f'Pix) < 0. (5.16) 

[atAh,kVhP{x) = at{x)AhA\x) 
+ / [al{x + ey)-al{x)]C{y)/^h,kVh{x + ey)dy, 



Next, 



JB 

where owing to Theorem II. 2^ for x E f^g/i+e the last term by magnitude is 
less than 

Ne / — ^ -dy, 

Jb P{x + ey) - 6/1 

which is less than Nep~^{x) if x G ^i2h.+2e since then 

p{x + ey) - p{x) > -e, p{x + ey) - {l/2)p{x) > {l/2)p{x) - e > 6h, 
p{x + ey)-6h> {l/2)p{x). 
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Hence in fti2h+2e 

al^h.kvi^ < [a1^h,kVhP+Nep~\ 
Then again owing to Theorem 11.21 in M 

ht5h,kvi^ < [bt6h,kVhP + Ne, -c'^v^^^ < -[c'^VhP + Ne, 

r < in^'^ + Ne. 

Coming back to (15.160 we conclude that in Qi2h+2e we have 

a^Af^^kvi'^ + btSh,kv'h^ - c^vl^^ + r{x) < Nep-\ (5.17) 

Now we are going to use Lemma 15.11 in order to replace A/j^/t and 5h^k 
in (|5.17p with D^ and -D;^, respectively. Of course, this time we take 
ri{x) = s~'^(^{x/e) and u = f /^ in Lemma 15. 1[ This lemma is stated only 
for vectors on the first basis axis. The reader understands that one can 
prove an appropriately modified statements for any vector I G B under the 
assumption that (j5.ip holds in the cylinder C2e+4:h,e,i centered at the origin 
with hight 2e + 4/i, base radius e, and axis parallel to I. Naturally, one 
can move such cylinders to be centered at any point xq. Observe that if 
3^0 £ ^2e+8h, then, for any / G A, all points of the cylinder xq + C2e+4h,£,i 
are at a distance not less than 

p{xo) -2e-2h>6h 

from Q'^. Hence, while applying Lemma |5. 11 to xq + C2£+ih,e,i we can take 
the constant N2 to be 

N{p{xo)-2£-8h)-\ 

where N is the constant from Theorem ll.2i Notice that 

p{xo) -2£-8h> l/2p{xo) 

if xo £ ^4e+i6h- Therefore, on r24£+i6/i by Lemma [5TT] we obtain 

\Di,vi'^ - 6h,kvP\ < Nhp-^ + Nh\e + h)e-\ 

\Df^vf - Ah,kVh^\ < Nhh^^p-^ + Nh\e + /i)e-^ 

These estimates and (j5.17|) yield (j5.15p . The lemma is proved. 

Assertion (i) of Lemma 15.21 that \v^ — Vh\ < N{e + h) in M can be 
improved for points at a fixed distance from $7^ if one uses the following 
lemma, in which ei, ..., e^ is the standard orthonormal basis in M . 

Lemma 5.3. Let u he a hounded Borel function on M'^. Assume that for all 
x,y £ [-{2h + e), 2h + e]'^ we have 

\u{x)-u{y)\ <Ni{\x-y\+h), 

\Ah,eM^)\<N2 i = l,...,d. 
Then for e > h 

\u^'\0)-u{0)\<N{£^ + h), (5.18) 
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where the constant N depends only on (^,d,Ni, and N2- 

Proof. We follow closely the argument in Remark 4.5 of [6]. For r > 
introduce 

Wr{x) = r^'^/'^ I Cr{x-y)u{y)dy, 
where 

Simple computations show that 

where (with no summation in i) 

Ci,r{x) = Ci{xr~^''^), Ci{x) = --l C{x - Xiei + sei)sds. 

Observe that since ^ is spherically symmetric, the support of Q is in B and, 
of course, d G Co°(K ). Also notice for the future that (no summation in i) 



\DfQA^)\ = ^}DtQ\{xr-'/') < l(|A'Cl + \DfC\){xr-y'). 
It follows that 



/¥* ^ rt* ^ 



* 1 



wtiO) - WriO) = J ^Dt[u*Q,smds. 

The support of Ci,s lies in B/^. Therefore, if < r < t < e^ and s S [r,t], 
we can use (|5.3p and obtain (no summation in i) 

= N^s'^/'WQU, + N,h\s'l^ + h)s''l^-^D%U,. 
Hence 

|u;j(0) - ^1-^(0)1 <N I {1 + h^s"^/"^ + h^s"'^) ds 

J r 

< N{t - r) + iV/iV"^/2 _^ AT/jS^-i^ 

where and below the constants N depend only on (^,d,Ni, and A^2- We 
combine this with 

\wriO) - um = I / CivMyr^^^) - n(0)] dy\ < Ni{r^/^ + h). 



Then we get 

\wt{0) - n(0)| < N{t - r) + 7V/iV~^/2 _^ jv/j3^-i _^ 7v(^i/2 _^ /j)^ 

which leads to (|5.18p if we take t = e^ and r = h^. The lemma is proved. 

Proof of Theorem 11.41 Take a constant n > 1 from Lemma 12.41 By 
Lemmas 15.2( 11) and 12.41 there is a constant A^ independent of a, h, and e 
such that for 

Ni := N{h + e + h^e'^) + Nh^e''^{e + h){e^ + h), 
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and sufficiently small h > we have 

<A>h + ^fc A,< - c"< + r < 

for all a G ^ in 17^, where k = (fih) V (4e+ 16/i). By the maximum principle 
V <wf^ + max(f — wf^)^ < vj^ + Ni^ + max(f — vjf )+ (5.19) 

in il.. By Theorem 11.31 we have |u — (7I < Nk in il. \ il^- Furthermore, by 
Lemma l5.2( i) we have \vj^ — Vh\ < N[e + h) everywhere which implies that 

l^i'^-^l <\vh-g\+N{e + h) 
and along with Theorem 11.21 yields that 

Ivl^^ -g\ <Nk + N{e + h) 

mQ\Qi^. Upon combining this with (I5.19P and observing that k < N{e + h) 
we get 

v<Vh + N{h + e + h^e-^) + Nh^e'^ie + h){e^ + h), 
which yields that in il for sufficiently small /i > we have 

v<Vh + Nh^/^ (5.20) 

if we set e = h"^'^. 
To prove that 

Vh<v + Nh^/^ (5.21) 

we reverse the roles of v^ and v. In i}^ we have 

KA> + btDi.v - c-v + r](^) < 0. (5.22) 

Furthermore, for x G il^ 

a'',Dlv^'\x)<[atDlv]^'^ + Ne sup |A>I, 

k,x+£B 

where by Theorem ll.3l the second term on the right is dominated by N{p{x) — 
e)~^, which is less than Np~^{x) if p{x) — e> (l/2)/9(x), that is if x G VL2e- 
Similarly one estimates h'^Di^v^^' , —d^v^^' and [Z"]^^-*. Then one concludes 
from (f5^ that 

a'^Df^'^ + hlDi^v^^^ - c°r;(^) + /" < Nep-^ (5.23) 

in ^2e- 

Next, as in ()5.6p for x G il2h+2£ we have 

\Dlv^'\x)-^h,kV^'\x)\<Nh' sup |A>(^)(y)|, 

where 

\Dlv^'\y)\ = \Dl[Dlv]^^\y)\ < Ne-\p{x) - h - e)-' < Ne-'p-\x) 
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if y G x + hB. Also as in (15.9 



where 



|A,^^^'^(a:) - 6h,kV^'>i^)\ < ^^^ sup |D;>^ 

x+hB 



\Dlv^'\y)\ < N{p{x) -h- e)~^ < Np-\x) 
if y € X + /li? and x G ri2/i+2e- Hence ()5.23p implies that 

in Vi2h+2e- At this point it is convenient to extend v outside Vl as g. By 
using Lemma [23] and the maximum principle as above we obtain that in W^ 

Vh < v^""^ + N{e + h'^e-^ + h) + sup{vh - v'^''^)+, 

where x = (/^^) V (2/i + 2e). Furthermore, in il. \ Q^ we have \vh — g| < Nx, 
which follows from Theorem II. 2( ii) and |u'^' — ffl < -^(x + s) ^ -^X; which 
follows from the fact that v G C^'^{Cl) by Theorem 11.31 Since x ^ iV(/i + e) 
we conclude that 

Vh <v^'^ + N{e + h^e-^ + h). 

The boundedness of the first derivatives of v implies that v^^' < v + Ne. 
Hence, 

Vh <v + N{e + h^e~^ + h), 

which yields ()5.2ip for e = h'^'^ and along with (I5.20p brings the proof of 
the theorem to an end. 
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